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ABSTRACT 

The axisymmetric 3-D MHD outflow of cold plasma from a magnetized and rotating 
astrophysical object is numerically simulated with the purpose of investigating the out- 
flow's magnetocentrifugal acceleration and eventual collimation. Gravity and thermal 
pressure are neglected while a split-monopole is used to describe the initial magnetic 
field configuration. It is found that the stationary final state depends critically on a 
single parameter a expressing the ratio of the corotating speed at the AlfVen distance 
to the initial flow speed along the initial monopole-like magnetic fieldlines. Several 
angular velocity laws have been used for relativistic and nonrelativistic outflows. The 
acceleration of the flow is most effective at the equatorial plane and the terminal flow 
speed depends linearly on a. Significant flow collimation is found in nonrelativistic 
efficient magnetic rotators corresponding to relatively large values of a >1 while very 
weak collimation occurs in inefficient magnetic rotators with smaller values of a < 1. 
Part of the flow around the rotation and magnetic axis is cylindrically collimated 
while the remaining part obtains radial asymptotics. The transverse radius of the jet 
is inversely proportional to a while the density in the jet grows linearly with a. For 
a >5 the magnitude of the flow in the jet remains below the fast MHD wave speed 
everywhere. In relativistic outflows, no collimation is found in the supersonic region 
for parameters typical for radio pulsars. All above results verify the main conclusions 
of general theoretical studies on the magnetic acceleration and collimation of outflows 
from magnetic rotators and extend previous numerical simulations to large stellar 
distances. 
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1 INTRODUCTION 

Plasma outflows from the environment of stellar or galactic 
objects, in the form of collimated jets is a widespread phe- 
nomenon in astrophysics. The most dramatic illustration of 
such highly collimated outflows may be perhaps found in rel- 
atively nearby regions of star formation; for example, in the 
Orion Nebula alone the Hubble Space Telescope (HST) has 
observed hundreds of aligned Herbig-Haro objects (O'Dell 
& Wen, 1994). In particular, recent HST observations show 
that several jets from young stars are highly collimated 
within 30 - 50 AU from the source star with jet widths 
of the order of tens of AU, although their initial opening 
angle is rather large, e.g., > 60°, Ray et al (1996). There 
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is also a long catalogue of jets associated with AGN and 
possibly supermassive black holes (Jones and Wehrle 1994, 
Biretta 1996). To a less extend, jets are also associated with 
older mass losing stars and planetary nebulae, (Livio 1997), 
symbiotic stars (Kafatos 1996), black hole X-ray transients 
(Mirabel & Rodriguez 1996), supersoft X-ray sources (Ka- 
habka & Trumper 1996), low- and high-mass X-ray binaries 
and cataclysmic variables (Shahbaz et al 1997). Even for 
the two spectacular rings seen with the HST in SN87A, it 
has been proposed that they may be inscribed by two pre- 
cessing jets from an object similar to SS433 on a hourglass- 
shaped cavity which has been created by nonuniform winds 
of the progenitor star (Burderi and King, 1995, Burrows et 
al 1995). 

In the theoretical front, the morphologies of collimated 
outflows have been studied, to a first approximation, in the 
framework of ideal stationary or time-dependent magneto- 
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hydrodynamics (MHD). In stationary studies, after the pio- 
neering 1-D (spherically symmetric) works of Parker (1963), 
Weber & Davis (1967) and Michel (1969), it was Suess 
(1972) and Nerney & Suess (1975) who first modelled the 
2-D (axisymmetric) interaction of magnetic fields with rota- 
tion in stellar winds, by a linearisation of the MHD equations 
in inverse Rossby numbers. Although their perturbation ex- 
pansion is not uniformly convergent but diverges at infinity, 
they found a poleward deflection of the streamlines of the so- 
lar wind caused by the toroidal magnetic field. Blandford & 
Payne (1982) subsequently demonstrated that astrophysical 
jets may be accelerated magnetocentrifugally from Keple- 
rian accretion disks, if the poloidal fieldlines are inclined by 
an angle of 60°, or less, to the disk midplane (but see also, 
Contopoulos & Lovelace (1994), Shu et al, 1994, Cao, 1997, 
Meier et al 1997). This study introduced the "bead on a rigid 
wire" picture, although these solutions are limited by the 
fact that they contain singularities along the system's axis 
and also terminate at finite heights above the disk. Saku- 
rai (1985) extended the Weber & Davis (1969) equatorial 
solution to all space around the star by iterating numeri- 
cally between the Bernoulli and transfield equations; thus, a 
polewards deflection of the poloidal fieldlines was found not 
only in an initially radial magnetic field geometry, but also 
in a split-monopole one appropriate to disk-winds, Sakurai 
(1987). The methodology of meridionally self-similar exact 
MHD solutions with a variable polytropic index was first 
introduced by Low & Tsinganos (1986) and Tsinganos & 
Low (1989) in an effort to model the heated axisymmetric 
solar wind. Heyvaerts & Norman (1989) have shown analyt- 
ically that the asymptotics of a particular fieldline in non 
isothermal polytropic outflows is parabolic if it does not en- 
close a net current to infinity; and, if a fieldline exists which 
does enclose a net current to infinity, then, somewhere in 
the flow there exists a cylindrically collimated core. Later, 
Bogovalov (1995) showed analytically that there always ex- 
ists a fieldline in the outflowing part of a rotating magne- 
tosphere which encloses a finite total poloidal current and 
therefore the asymptotics of the outflow always contains a 
cylindrically collimated core. In that connection, it has been 
shown in Bogovalov (1992) that the poloidal fieldlines are de- 
flected towards the polar axis for the split monopole geome- 
try and relativistic or nonrelativistic speeds of the outflowing 
plasma. Sauty & Tsinganos (1994) have self-consistently de- 
termined the shape of the fieldlines from the base of the out- 
flow to infinity for nonpolytropic cases and provided a sim- 
ple criterion for the transition of their asymptotical shape 
from conical (in inefficient magnetic rotators) to cylindrical 
(in efficient magnetic rotators). They have also conjectured 
that as a young star spins down loosing angular momentum, 
its collimated jet-type outflow becomes gradually a conically 
expanding wind. Nevertheless, the degree of the collimation 
of the solar wind at large heliocentric distances remains still 
observationally unconfirmed, since spacecraft observations 
still offer ambiguous evidence on this question. Another in- 
teresting property of collimated outflows has emerged from 
studies of various self-similar solutions, namely that in a 
large portion of them cylindrical collimation is obtained only 
after some oscillations of decaying amplitude in the jet-width 
appear (Vlahakis & Tsinganos 1997). Radially self-similar 
models with cylindrical asymptotics for self-collimated and 
magnetically dominated outflows from accretion disks have 



been constructed in Ostriker (1997). All existing cases of 
self-similar, jet- or, wind-type exact MHD solutions can be 
unified by a systematic analytical treatment wherein the 
available today examples of exact solutions emerge as special 
cases of a general formulation while at the same time new 
families with various asymptotical shapes, with (or without) 
oscillatory behaviour emerge as a byproduct of this system- 
atic method (Vlahakis & Tsinganos, 1998). Altogether, some 
general trends on the behaviour of stationary, analytic, ax- 
isymmetric MHD solutions for MHD outflows seem to be 
well at hand. 

However, observations seem to indicate that jets may 
inherently be variable. Thus, time-dependent simulations 
may be useful for a detailed comparison with the observa- 
tions. Uchida & Shibata (1985) were the first to perform 
time-dependent simulations and demonstrate that a verti- 
cal disk magnetic field if twisted by the rotation of the disk 
can launch bipolar plasma ejections through the torsional 
Alfven waves it generates. However, this mechanism applies 
to fully episodic plasma ejections and no final stationary 
state is reached to be compared with stationary studies. 
Similar numerical simulations of episodic outflows from Ke- 
plerian disks driven by torsional Alfven waves on an initially 
vertical magnetic field have been presented by Ouyed & Pu- 
dritz (1997a,b). Goodson et al (1997) have proposed a time- 
dependent jet launching and collimating mechanism which 
produces a two-component outflow: hot, well collimated jet 
around the rotation axis and a cool but slower disk-wind. 
Stationary MHD jet-type outflows have been found in the 
study of Romanova et al (1997), as the asymptotic state 
of numerical simulations wherein a small initial velocity is 
given in the plasma in a tapered monopole-like magnetic 
field. Numerical viscosity however results in nonparallel flow 
and magnetic fields in the poloidal plane in the limited grid 
space of integration. Washimi & Shibata (1993) modelled 
axisymmetric thermo-centrifugal winds with a dipole mag- 
netic flux distribution B 2 (#) oc (3cos 2 # + 1) on the stel- 
lar surface (and a radial field in Washimi 1990). In this 
case the magnetic pressure distribution varies approximately 
as oc Bp sin 2 9 such that it has a maximum at about 
cos~ 1 20 o « —1/3, or, do « 55°. As a result, the flow and 
flux is directed towards the pole and the equator from the 
midlatitudes around 9 . The study was performed for uni- 
form in latitude rotation rates and up to 60 solar radii in 
the equatorial plane. Bogovalov (1996, 1997) modelled nu- 
merically the effects of the Lorentz force in accelerating and 
collimating a cold plasma with an initially monopole-type 
magnetic field, in a region limited also by computer time, 
i.e., the near zone to the central spherical object. 

This paper presents an extension of the previous results 
to large distances from the star by using a new method for 
the continuation of the small simulation box solution to very 
large distances. It also examines the efficiency of magnetic 
rotators of various strengths in transforming rotational en- 
ergy to directed kinetic energy and how a wide range of rota- 
tion rates affects the poloidal geometry of a magnetic field. 
In order to achieve these objectives, the rather complicated 
nature of the problem requires that we start by limiting the 
investigation to the simplest model of cold plasma flow in 
a uniform in latitude monopole magnetic field along which 
there is an outflow with velocity V , as the initial condition. 
This simplification allows us to better study and understand 
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the nonlinear effects of the magnetocentrifugal forces alone 
in shaping the final stationary configuration. At the same 
time, however, it does not allow us to perform a direct com- 
parison of the obtained results with observations. Clearly 
this is not the goal of the present paper since for such a 
comparison we need to include gravity and thermal pres- 
sure in our computations. Such a study has been already 
performed in the context of the solar wind and it will be 
presented elsewhere. 

The paper is organised as follows. In Sees. 2 and 3, 
the initial configuration used together with the method for 
the numerical simulation in the nearest zone is discussed. In 
Sees. 4 and 5 the analytical method for extending the inte- 
gration to unlimited large distances outside the near zone 
are briefly described. In Sees. 6 and 7 we discuss the results 
in the near zone containing the critical surfaces and in the 
asymptotic regime of the collimated outflow, for a uniform 
rotation. In Sec. 8 a rotation law appropriate for an accre- 
tion disk is used, while in Sec. 9 we briefly discuss results of 
a relativistic modelling. A brief summary is finally given in 
the last Section 9. 



obtained from a monopole-like solution by a simple revers- 
ing of the magnetic field in one of the hemispheres. 

(2) In considering the plasma flow at such large distances 
it is natural to neglect gravity. Thermal pressure however 
may play an important role even at large distances from the 
central object (Bogovalov 1995). But in this paper we are 
interested to isolate the effects arising purely from the mag- 
netic field. Thus, to make our analysis as simple as possible, 
in this paper we neglect thermal pressure too. 

(3) We are also interested to study the magnetocentrifu- 
gal acceleration of the plasma. This acceleration process can 
be studied with a monopole like magnetic field regarded as 
a first approximation to the more realistic acceleration in a 
dipole-like magnetic field with open field lines. Gravity and 
thermal pressure can be neglected in this case (Michel 1969). 

The previous claims (i) - (iii) are the result of a rather gen- 
eral analysis. Among the main goals of the present work is 
to verify these general conclusions in the context of the split 
monopole model and assumptions (1) - (3). 



2 THE MODEL OF A ROTATOR WITH A 

MONOPOLE-LIKE MAGNETIC FIELD AND 
THE OBJECTIVES OF THIS PAPER 

A general analysis of the asymptotical properties of nonrel- 
ativistic or relativistic magnetized winds has been already 
performed, e.g., Heyvaerts and Norman (1989), Chiueh et al 
(1991), Bogovalov (1995). The main conclusions from such 
studies can be summarized as follows: 

(i) At large distances from the central source the poloidal 
magnetic field is similar (although not exactly the same) to 
a split-monopole field, 

(ii) there exist cylindrically collimated and radially ex- 
panding field lines, while the total electric current enclosed 
by any magnetic surface is nonzero, 

(iii) several physical quantities accross the jet can be ex- 
pressed by simple formulas, under certain conditions. 



The model of an axisymmetric rotator with a monopole- 
like magnetic field was first used by Michel (1969) for the 
investigation of the cold plasma flow in the absense of grav- 
ity in a prescribed poloidal magnetic field. Later this model 
was used by Sakurai (1985) in an attempt to solve selfcon- 
sistently the problem of the nonrelativistic plasma outflow 
from a stellar object. This model may be used to the study 
of plasma outflow from the magnetospheres of various cos- 
mic objects under the following conditions: 

(1) We are interested in the plasma flow at large distances 
where all magnetic field lines are open. In this case, the field 
of any axisymmetric rotator, no matter what is the nature 
of the central object, becomes the field of the so-called split- 
monopole in the far zone. In other words, the condition that 
should be fulfilled is that the distance to the central source 
should be much larger than all the dimensions of the cen- 
tral source (see also Heyvaerts & Norman 1989, Bogovalov 
1995). The solution for the split monopole field can be easily 



3 THE PROBLEM IN THE NEAREST ZONE. 

To obtain a stationary solution of the problem in the nearest 
zone of the star containing the critical surfaces, it is needed 
to solve the complete system of the time-dependent MHD 
equations and look for an asymptotic stationary state. In 
order to isolate the effects of the magnetic field in determin- 
ing the shape of the streamlines, we shall neglect gravity 
and thermal pressure gradients, as discussed above. With 
these simplifications, the flow of the nonrelativistic plasma 
is described by the set of the familiar MHD equations, 
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dt dr dz 8irpr 2 dr 



9 irB v f + 



Yl + j^_ ( 9B r 

r 4np \ dz 



dB, 
dr 



(7) 



where we have used cylindrical coordinates (z, r,ip), p is the 
density, V the flow field and B the magnetic field with a 
poloidal magnetic flux denoted by ip(z,r). 

A correct solution of the problem requires a specifica- 
tion of the appropriate boundary conditions at some spher- 
ical boundary R = R of the integration. 

1. A constant plasma density p at R = R . 

2. A constant total plasma speed Vo in the corotating frame 
of reference at R = R a , V ( 2 r _ o) + V ( 2 Z}0) + {V (f , o) - fir ) 2 = V 2 

3. A constant and uniform in latitude distribution of the 
magnetic flux function ip — tp at R = R . 

4. Finally, the continuity of the tangential component of the 
electric field across the stellar surface in the corotating frame 
gives the last condition, (V (¥!>0 ) - Q,r )B {Pt0) - V^ p>a )B^^ = 
0. 

We shall use dimensionless variables, Z = z/R a , X = 
r/R a , t = tVo/Ra where R a is the Alfven spherical radius 
of an initially radial, monopole-like, nonrotating magnetic 
field and define the dimensionless parameter 

■ (8) 
This parameter a characterizes the influence of the magnetic 
field and rotation on the acceleration and collimation of the 
plasma. It is proportional to the time the plasma spends in 
the subAlfvenic region in each period of rotation. Note that 
although the governing equations (1-7) do not depend on a, 
the final solution does depend on a through the boundary 
conditions at the base of the integration R = R (condition 
4). 

As we shall see, the solution in the nearest zone will 
relax after sufficient time to a stationary state. This station- 
ary state will be next used as the input for specifying the 
boundary conditions in the super-fast magnetosonic region. 
In this way we shall be able to obtain a complete solution of 
the stationary problem, from the base up to large distances 
downstream. 



4 THE STATIONARY PROBLEM 

Below we shall consider that the plasma may be relativis- 
tic, or nonrelativistic and as before we shall neglect gravity 
and thermal pressure. By U p = 'yVp/c and U v = jVtp/c, we 
denote the poloidal and azimuthal 4-speeds where V v is the 
azimuthal and V p the poloidal components of the velocity 
while 7 is the Lorentz factor of the plasma. In the follow- 
ing subsection we review the basic quantities which remain 
invariant along a poloidal streamline ip = const, and ex- 
press momentum balance along such a poloidal streamline. 
In the next subsection we adopt a new coordinate system for 
dealing with the transfield equation expressing momentum 
balance across the poloidal streamlines. 



4.1 MHD Integrals 

As is well known, the stationary MHD equations admit four 
integrals. They are: 



(a) The ratio of the poloidal magnetic and mass fluxes, 
cF(V>) 

(/3) The total angular momentum per unit mass L{ijj), 

rcU v - cFrB v = L(ip) . (10) 

(7) The corotation frequency Q(ip) in the frozen-in condi- 
tion 



cU v B p — cUpBy, = rjB p Q(ip) ■ 



(11) 



(S) The total energy c 2 W(ip) in the equation for total en- 
ergy conservation, 

7c 2 - cF(i>)rSl{il))B v = c 2 W(tP) . (12) 



4.2 The transfield equation in the coordinates 

Momentum balance across the poloidal fieldlines is expressed 
by the transfield equation which determines and their shape. 
This is a rather complicated nonlinear partial differential 
equation of mixed elliptic/hyperbolic type. For analysing 
the behaviour of the plasma at large distances, it occured 
to us that it is rather convenient to work with this trans- 
field equation in an orthogonal curvilinear coordinate system 
(tf>, rj) formed by the tangent to the poloidal magnetic field 
line fj = p and the first normal towards the center of curva- 
ture of the poloidal lines, rp = V^/|V^| (Sakurai, 1990). A 
geometrical interval in these coordinates can be expressed 



(dr) 2 = g^dtp 2 + g^drj 2 + r 2 dip 2 . 



(13) 



where g^,g v are the corresponding line elements or compo- 
nents of the metric tensor. 

If T !J is the energy-momentum tensor of the plasma 
flow (pV) and electromagnetic field (E, B) (Landau & Lif- 
shitz 1975), the equation dT^ /dx k = (with covariant 
derivatives) has the following form in the coordinates (tp, rj), 



d 


'B 2 -E 2 ' 


1 dr 


dtp 


8tt 


r dip 



Bt - E 1 



4tt 



1 dg v 

9v dip 



pVp 2 



B 2 -E 2 
4% 



0. 



(14) 



The first term in this equation is the gradient of the 
pressure of the electromagnetic field, while the second is the 
sum of the inertial terms due to the motion of the plasma 
in the azimuthal direction and also due to the tension of the 
toroidal magnetic field. To better understand the physical 
meaning of the last term, note that since 77 is perpendicular 
to tp we have, 



(15) 



1 dg v = 1 
g v dip rRcBp ' 

where R c is the radius of curvature of the poloidal magnetic 
field lines, with R c positive if the center of curvature is in 
the domain between the line and the axis of rotation and 
negative in the opposite case. 

With this expression, Eq. (|14|) becomes, 
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[r\Bl-E 2 )] 



1 dr 



U 2 B P 



4nr dtp U p F(ip) 



p p -FW(l-(#) a R] 
4-KrR c F(ip) 



0. (16) 



From this equation, it may be seen that the last term is the 
sum of the inertia of plasma and fields, connected with the 
motion of the plasma and Poynting flux along the poloidal 
field line and tension of the poloidal field line. 



4.3 The transfield equation for nonrelativistic 
plasmas with gravity and thermal pressure 
included 

For completeness of the picture we present briefly here 
the transfield equation for a nonrelativistic plasma flow to 
demonstrate that our method of solution of the stationary 
problem in the hyperbolic region can be applied directly 
to this case too. The energy momentum conservation equa- 
tion in the presence of gravity in the nonrelativistic limit is 
modified as dT^ h jdx h = -pd^/dx^, where T^ k is again 
the energy-momentum tensor, $ = —GM/R is the gravi- 
tational potential of the star with mass M and G is the 
gravitational constant. This equation has the following form 
in our curvilinear coordinates if some thermal pressure P is 
also included, 



d_ 

dip 



P + 



8tt 



1 dr 

r dip 



B v 

47T 



1 dg v 

9v dip 



= -P 



9£ 

dip 



(17) 



All other equations, describing the flow of plasma along 
field lines will be the same except the energy conservation 
equation which is modified as follows 

^- + T ^—- + ^~cF{iP)rQ(iP)B^ = W(iP)c 2 . (18) 
Z o — 1 p 

where <5 is the polytropic index of the plasma. 

5 THE SOLUTION IN THE FAR ZONE. 

It is convenient to solve the transfield equation in the sys- 
tem of the curvilinear coordinates introduced above. The 
unknown variables are z(jj,ip) and r(rj,ip). Therefore we 
need to know the quantities g v , g^,, r^, z^,, r v , z v , where 
r v = dr/drj, z v — dz/drj, = dr/dip, z^ = dz/dip. 

First, the metric coefficient g v is obtained from the 
transfield equation (H), 



g v = exp ( / G(r], ip)dip) . 



where 
G(ri,ip) 



(19) 



g 

Oil' 



[^(B 2 -E 2 )] 



1 dr_ 

r dip 



K 2 -^-£ 2 )] 



[PV P 2 - M B p - E2 )] 



(20) 



for a cold plasma. For a nonrelativistic flow with finite ther- 
mal pressure and gravity the function G will have the form 

G(v,ip) = 



(21) 



gg [^b 2 +£\+^z!m ^ ~ Egg] 

[PV P 2 - ^Bl] 

The lower limit of the integration in Eq. (jl9|) is chosen to 
be such that the coordinate rj is uniquely defined. In this 
way rj coincides with the coordinate z where the surface of 
constant r\ crosses the axis of rotation. 

Second, the metric coefficient g^ is given in terms of the 
magnitude of the poloidal magnetic field by, 



9i> = 



rBv 



(22) 



To obtain the expressions of r,p, r v , z v we may use 
the orthogonality condition 



(23) 



and also the fact that they are related to the metric coeffi- 
cients g v and as follows, 



2 

9n 



2 , 2 
' \) I" Zjj , 



2 2 . 2 



(24) 
(25) 



Thus, by combining the condition of orthogonality ( |23| 
and equations (^i|) and (^B[) the remaining values of r v , z, 
are obtained, 



9i> 
T^9n 



(26) 



(27) 



with g v calculated by the expression (h9|) . For the numerical 
solution of the system of equations (fjq - a two step Lax- 
Wendroff method is used on a lattice with a dimension equal 
to 1000. 

Equations (^ - ^?j) should be supplemented by appro- 
priate boundary conditions on some initial surface of con- 
stant rj. The equations for r^ , z^ defining this initial surface 
in cylindrical coordinates are as follows 

dr - B > 



rBl 



dz 
£hp 



Br 

VbI 



(29) 



We need to specify on this surface the integrals 
F(ip), L{ip), Q(ip) and W(ip) as the boundary conditions for 
the initial value problem. To specify the initial surface of 
constant rj and the above integrals, we use the results of the 
solution of the problem in the nearest zone when a station- 
ary solution is obtained for the time-dependent problem. 



6 RESULTS IN THE NEAREST ZONE FOR 
UNIFORM ROTATION, Q.(tp) = flo 

As the star starts rotating, in the initially radial magne- 
tosphere an MHD wave propagates outwards carrying the 
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Cold plasma 
a = 2. 4780184 





i 1 / y 1 / 1 / 
/ / / / / 


mfffi 
\\\\\\ / 




J ' 1 \ / 
j i | 

J ' \ 
nil i j 

[ l / 




f / 
' 




I i j j 
F'l' I'! I 1 ''' 




ffllf/'/// / 


iJ/////^' 


Wiit/7 

l 


-rr, n , i 



2 4 6 



Figure 1. A near zone snapshot on the poloidal plane (Z, r) from 
the simulation showing the change of the shape of the poloidal 
magnetic fieldlines from an initially uniform with latitude radial 
monopole and before a stationary state is reached. Distances are 
given in units of the Alfven radius R a and time in units of the 
Alfven crossing time t a = R a /V a 



effect of the rotation and deflecting the fieldlines polewards 
by the Lorentz force (Fig. 1). 

At times sufficiently long after the shock wave has 
reached the end boundary of the simulation (t » 1), a 
final equilibrium state is reached (Figs. 3). In this station- 
ary state, the poloidal magnetic field and the plasma density 
are increased along the axis because of the focusing of the 
field lines towards the pole. 

In Figs. (3) the distances are given in units of the ini- 
tial Alfvenic radius R a . Thick lines indicate Alfven and fast 
critical surfaces while thin lines the poloidal magnetic field. 
At t — the Alfven surface is spherical at R = 1. Dotted 
lines indicate poloidal currents. 

With the velocity maintaining the constant initial value 
V along the Z-axis, the ratio B p /p remains constant along 
this axis (tp — 0). As a result, the Alfven speed at a 
given point of the Z-axis, Va = Bp/yJAitp, increases as the 
square root of the density increases there because of focus- 
ing. It follows that the Alfven transition at Ra(0 — 0) oc- 
curs further downstream where V (Ra) = Va{Ra), i.e., at 
Ra(0 — 0) > 1. As we move meridionally off the polar axis 
toward the equator on the other hand, the bulk flow speed 
is increased because of the magnetocentrifugal acceleration; 
it thus hits the Alfven value earlier than it does on the axis, 
i.e., V(6) = V A (9), at R A {0) < Ra(9 = 0). Finally, on the 
equator the flow speed increases rapidly with the result that 
the flow becomes much earlier super Alfvenic there in com- 
parison to the polar axis. The degree of elongation of the 
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Figure 2. Maximum height Z max of fast critical surface as a 
function of a, in units of the base radius R a . 

Alfven surfaces along the symmetry axis increases with the 
value of q. 

The behaviour of the fast surface has similarities and 
differences with the shape of the Alfven surface. First, both 
these critical surfaces coincide at 6 — 0. As we move merid- 
ionally off the polar axis toward the equator on the other 
hand, their shape becomes different. This is due to the con- 
tribution of the azimuthal field, since in this case of cold 
plasma Vf = (Bp + B^/Aivp. Thus, for small colatitudes 
8, Vf increases due to the contribution of B<f, while also V 
slightly increases because of the magnetocentrifugal acceler- 
ation. However, initially Vf increases faster and therefore the 
fast transition is postponed further downstream. At larger 9 
however, the competing flow speed V increases faster than 
Vf does and the fast transition occurs closer and closer to 
the origin. 

The dependance of the maximum height Z max of the 
fast critical surface as a function of a is interesting too (Fig. 
2). As a increases, Z max increases rapidly. That means that 
for sufficiently fast rotators with a ^ 5, Z max goes to infin- 
ity and the flow may stay sub-fast all the way to large dis- 
tances. This result may have some important implications 
in the general theory of MHD winds. It clearly indicates 
that a super-fast stationary solution may not be obtained 
for all sets of parameters. For large a no super-fast station- 
ary state is found and thus the equilibrium is vulnerable to 
instabilities. Such a situation may take place, for example, 
in young rapidly rotating stars. Outflows from classical ac- 
cretion disks (Shakura & Sunyaev 1973) are also expected 
to have such large values of a, since the Alfven and sound 
speeds in such disks are much smaller than their respective 
Keplerian speed. Thus, since V a ~ V s << Vk and r a > tk, 
we have 

Q,r a Qr K , 
a = ~V7 ~~Vk ~ ' 
It follows that outflows from rapidly rotating stars and 
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Figure 3. Sequence of shapes of the poloidal field lines by increasing the magnetic rotator parameter a from a = 0.5 (solar wind 
type slow magnetic rotator) to a = 4.5 (fast magnetic rotator). The initial nonrotating monopole magnetic field has a spherical 
Alfven surface located at R = 1 and all distances are given in units of the Alfven radius R a with the base located at R = 0.5. 
Dotted lines indicate poloidal currents. Thick lines indicate Alfven and fast critical surfaces. At t = the Alfven surface is spherical 
at R = 1. 
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thin classical accretion disks may rather produce jets with 
characteristics which qualitatevely differ strongly from those 
of laminar jets usually discussed in the literature. This result 
is certainly obtained under the simplifying assumptions of 
the present study, and further investigations are necessary 
to explore the possibilities which may arise in this regime of 
outflow. Up to now all numerical simulations avoided this 
problem by artificially taking a ~ f in order to obtain a 
stationary solution (Romanova et al 1997), or, no stationary 
solution was obtained at all (Ouyed & Pudritz 1997, see 
Fig. 8 in this work). Fig. (3) then may indicate why no 
simulation has succeeded to produce stationary supersonic 
jet from classical accretion disks up to now (Ferreira, 1997) . 



7 RESULTS IN THE FAR ZONE 

The characteristics of magnetized outflows at large distances 
from the central object have been studied by Heyvaerts & 
Norman (1989), Chiueh et. al. (1991) and Bogovalov (1995). 
These studies have concluded that a stationary axisymmet- 
rically rotating object ejecting magnetized plasma always 
produces a jet collimated exactly along the axis of rotation, 
if the following conditions are satisfied: 

1. The flow is nondissipative. 

2. The angular velocity of rotation is nonzero everywhere 
(actually, if Q — on some field lines, the same conclusion 
remains valid). 

3. The total magnetic flux reaching infinity in any hemi- 
sphere of the outflow is finite. 

4. The polytropic index of the plasma 5 > 1. 

In such an outflow the density of the poloidal electric current 
is nonzero in the region of the collimated flow. Conversely, 
in the region of noncollimated field lines the density of the 
poloidal electric currents equals to zero. This condition can 
be expressed by the constancy of the quantity r 2 Q(ip)B p /U p , 
which we shall call the Heyvaerts-Norman integral. 



7.1 Efficiency of the magnetic rotator 

In the Weber & Davis (1967) model of a magnetized equato- 
rial wind the terminal equatorial speed is superfast with the 
fast mode surface placed at some finate distance from the 
star. If the initial velocity of a cold plasma is equal to zero 
on the surface of the star, the fast critical point is at infinity 
and we obtain the so-called Michel's (1969) minimum en- 
ergy solution. This solution is valid when the monopole-like 
poloidal magnetic field is slightly disturbed by the plasma. 
In this case the asymptotic speed at infinity on the equator 
is Voo = (3/2) V a , while the total specific angular momentum 
in the system in units of R a V a is, L = Qr^/R a V - We recall 
that R a is the radius of the initial spherical Alfven surface 
while r a is the Alfvenic radius on the equator at a given 
angular velocity. An important physical quantity in mag- 
netized outflows is the magnetic rotator energy, Emr, the 
product of the total specific angular momentum and £1. The 
basal Poynting energy defined as the ratio of the Poynting 
flux density S z per unit of mass flux density pV z is approxi- 
mately equal to -Emr i/at the base of the outflow the radius 
of the jet is much smaller than the Alfven radius and also 
the Alfven number there is negligibly small. 



es 4.00 




3.00 



Figure 4. For a magnetic rotator we plot the rotational losses 
ctL per particle, as a function of a (solid line). For comparison, 
the corresponding rotational losses for Michel's minimum energy 
solution are also plotted (dashed line). 



The energy losses of the magnetic rotator per particle 
in units of V 2 are 



-Emr 



V? 



aL . 



(30) 



This formula is widely used in the theory of the rotational 
evolution of stars. But the parameter a goes to oo in Michel's 
minimum energy solution. The plasma is strongly collimated 
under this condition and therefore Michel's solution is not 
valid. In this case it is important to know how strongly the 
effect of collimation affects those frequently used Michel's 
energy losses. The energy losses aL per particle on the equa- 
tor for Michel's monopole-like solution and our calculations 
are compared in Fig. kl 

The decceleration rate in our solution is less than in 
Michel's solution. This is due to the collimation of the 
plasma. The poloidal magnetic field near the equator de- 
creases and it leads to a reduction of the decceleration rate. 

It is interesting that in spite of the decrease of the dec- 
celeration rate the terminal velocity of the plasma near the 
equator increases in our solution. The terminal speed Voo 
as a function of a is plotted in Fig. (^). The terminal ve- 
locity is calculated on the equator at the dimensionless dis- 
tance X=500. For the case of a fast magnetic rotator (Michel 
1969), the dependance of Voo on a is, 



Vo 



1 f ^RjBg .1/3 

Vo M 



v 2 / 3 



i.e., it goes like a 2//3 . This increase of the terminal velocity 
on the equator in comparison to that in Michel's solution is 
due to the collimation of the plasma. In the ollimated flow 
there exist strong gradients of the toroidal magnetic field 
which additionally accelerate plasma. This means that in 
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Figure 5. The terminal velocity Voo/Vo as a function of a 
(solid line). For comparison, the corresponding terminal speed in 
Michel's minimum energy solution is also plotted (dashed line). 
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Figure 6. Efficiency e of the magnetic rotator as a function of 
a (solid line). For comparison, Michel's solution has e = 1/3 
(dotted line). 



collimated ouflows we have a more effective acceleration of 
the plasma by the magnetic rotator. 

The efficiency e of the magnetic rotator in transforming 
part of the base Poynting flux to poloidal kinetic energy at 
infinity is measured as the difference of the poloidal kinetic 
energies at infinity and at the base normalized to the total 
energy E (in units of V D 2 ). The poloidal kinetic energy at 
infinity in Michel's solution in units of V 2 , E%£ 1 is 

t _E_ Q*rl 
3 3V? ' 

or, 




The efficiency e is plotted in Fig. (g) for our solution and 
(solid) and for Michel's solution (dashed). We see that due 
to collimation, the magnetic rotator becomes a very effective 
machine for plasma acceleration. 

7.2 The radius of the jet 

The dependence of the magnitude of the poloidal magnetic 
field and density on the cylindrical distance r becomes espe- 
cially simple if we assume for convenience that the following 
additional conditions are met in the jet, namely that the 
integrals W (?/>), L(ijj), F(tp), Q(4>) an d the terminal velocity 
Vj are constants and do not depend on ip and also that 
Vj 3> V a (0), where V a (0) is the Alfvenic velocity on the axis 
of rotation. In such a case, Eqs. ( |l6| ) - give an approxi- 
mate estimate of the dependence of the magnetic field on r 
(Bogovalov 1995), 

B P (r) _ p(r) _ 1 

Bp(0) P (0) (1 + (r/R^ry ( ' 

where Rj is the radius of the core of the jet 



Rj -T 1+ KW ) ~n> (32) 

with C s (0) the sound velocity along the jet's axis and B p (0), 
p(0) the magnetic field and the density of the plasma on the 
axis of the jet, respectively. It is evident that the poloidal 
magnetic field and the density remain practically constant 
up to distances of order Rj and then decay fast like 1/r 2 
outside the jet's core. 

Among the main goals of the present paper is the veri- 
fication of these conclusions, together with the investigation 
of the process of jet formation and the development of meth- 
ods for the calculation of the characteristics of the plasma 
in the jet. 

In Fig. (Q) the poloidal and azimuthal components of 
the magnetic field are plotted together with the magnetic 
flux enclosed by a cylindrical distance X. The poloidal mag- 
netic field B P (X)/ B (solid line) is given in units of its ref- 
erence value B , corresponding to the magnetic field at the 
symmetry axis r — and some reference height Z a — 60825 
for a = 0.5 and Z = 5000 for a — 3. The asymptotic 
regime of the jet is achived at these distances. The intensity 
of the poloidal magnetic field drops dramatically by more 
than 50% with respect to its value at the axis X = within 
a distance of the order of the collimation radius Xj = 1/a. 
The dashed curve gives the analytically predicted solution 
for the poloidal mag netic field B P (X)/B , Eq. (|l|). Appar- 
ently, for small values of a there is a good agreement be- 
tween the calculated and analytically predicted values of 
the poloidal magnetic field. For a — 3 the agreement be- 
tween the analytical prediction and numerical calculations 
is worse than for a — 0.5. This is natural because the con- 
dition Vj 3> V a (0) under which the analytical prediction is 
valid becomes not so good fulfilled for a — 3 as for a = 0.5. 

The strength of the azimuthal magnetic field B<f, (X) / B 
(dots) obtains a maximum value at the radius Xj. For 
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Figure 7. Variation with the dimensionless cylindrical distance X of the enclosed magnetic flux ip(X) (dot-dashes), strength of the 
azimuthal magnetic field B l f > (X)/ B (dots) and poloidal magnetic field B P (X) / B (solid) for a slow magnetic rotator a = 0.5 and a 
faster magnetic rotator a = 3. With dashes the analytically predicted solution for the poloidal magnetic field B p (r)/B is also shown. 



X < Xj we have approximately conditions similar to those 
corresponding to a uniform current density wire and there- 
fore B$ oc X. On the other hand, for X > Xj, conditions 
like those existing outside of a uniform current density wire 
exist and therefore Bj, oc X~ . 



7.3 The Heyvaerts & Norman enclosed current 

The enclosed electric current by a poloidal field line ip = 
const, is plotted in Fig. (|^a) for a = 3. We may distinguish 
three regimes. Close to the axis, < ip < 0.1 the current in- 
creases with ip since there we have conditions similar to those 
corresponding to some uniform current density wire. A new 
regime appears when we are outside the collimated region 
where the enclosed current reaches a plateau, 0.1 < ip < 0.5. 
We shall call this regime the Heyvaerts-Norman regime, 
since there we have conditions similar to those existing out- 
side a uniform current density wire (XB,f, = const.), a sit- 
uation described by Heyvaerts & Norman (1989). Finally, a 
third domain exists in 0.5 < ip < 1 where again the enclosed 
current increases. It seems that this contradicts the expected 
behaviour, but as we show below it is due to very slow de- 
crease of this part of the current and we do not achieve final 
asymptotics. 



7.4 Logarithmic collimation asymptotically 

One may use the coordinates (n, ip)to show the dis- 
tribution of the electric currents, Figs. (pis). In this space, 
ip = corresponds to the symmetry axis X = 0, ip = 1 to 
the equator and rj = approximately to the source surface, 
while a poloidal fieldline corresponds to some vertical line 
ip = const. Fig (^b) demonstrates that there is an electric 
current near the axis. Then the Heyvaerts & Norman region 
of constant XB V is formed. And finally the region where 
XB V increases is observed. It is important that the excess 
of the electric current in this last region decreases with rj. 
But this decrease occurs in a logarithmic scale. The solution 



goes to its asymptotic form, but this is done very slowly, in 
a logarithmic scale. This behavior can be demonstrated also 
by a simple analysis of the transfield equation. 

In a supersonic flow, the motion of every parcel of 
plasma is controlled by the initial conditions and the forces 
affecting the plasma. Here we are interested in the collima- 
tion process in the region of radially expanding field lines. 
In this region forces due to the poloidal magnetic field and 
the inertia of the plasma due to motion in the azimuthal 
direction can be neglected in the transfield equation in com- 
parison to forces arising from the toroidal magnetic field. 
The terms Bp/8ir and pV£ can be neglected since they drop 
with distance r as 1/r 4 , while the terms corresponding to 
the azimuthal magnetic field B^ drop as 1/r 2 . Then, the 
transfield equation Eq. (XL&) is simplified as follows, 



9 (rB v f 



dip 



ATvrR c F(tp) 



0. 



(33) 



It may be seen from this equation that the curvature 
radius of a field line of the poloidal magnetic field is de- 
fined by the tension of the toroidal magnetic field. Let us 
estimate how the curvature radius of the poloidal magnetic 
field changes with distance assuming for simplisity that the 
poloidal magnetic field expands radially and hardly depends 
on the polar angle 8 near equator. At large distances, the 
frozen in condition gives for the electric current 



rBu 



r 2 flB p 
V v 



(34) 



Inserting this in (g3|) and assuming that the dependence of 
Bp on 8 is weak, we get 



Vv 



v P 2 



rR c Fc' 



(35) 



The curvature radius of the field line is defned by the 
equation 



dl 



(36) 
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Figure 9. Shape of the poloidal magnetic field lines for a magnetic rotator parameter a = 3 plotted up to a radius 5000 times the initial 
base radius. In (a) the poloidal field lines are plotted in a logarithmic scale, and give the erroneous impression that all fieldlines are 
focused towards the system's symmetry axis. In the linear scale of Fig. (b) however, a cylindrical jet is formed around the symmetry 
axis while all other field lines go to straight asymptotes which fill all space. 



where dl is the element of the length of the fied line. Since 
to a first approximation dl — dR, where R is the spherical 
distance to the center of the rotating object, we get for the 
equation of a field line, 



dO 
dR 



cFBpR 2 D. 2 . 
RVi Sm26 - 



(37) 



For a radially expanding magnetic field B P R 2 is constant 
along a magnetic field line. Therefore, the turn angle of the 
field line A9 depends logarithmically on R, 

Vp Ro 
where _Ro is the initial distance. According to (|34|) the change 
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of the electric current enclosed by a field line is connected 
with the divergence of the poloidal magnetic field from the 
purely radial shape. To have variation of rB^ with ip it is 
necessary to turn some field lines at some angle. Such a 
turn occurs at the exponentially large distances defined by 
the expression 



R — Ro exp ( 



Vp 3 A6» 



cFB p (RQ) 2 sin 26 



) > Ro exp 



A6»tan0 



(39) 



because for the super fast magnetosonic plasma we have 



V p > cFB 2 v /B p , 



(40) 



which taking into account ( |34j ) can be rewritten as V p > 
cFB p r 2 Q. 2 . In other words, the magnitude of the exponent 
in Eq. ( pj9| ) is larger than A6 tan 6/2, a large number near 
the equator. Apparently, two points at small angular dis- 
tance A6 on the same fieldline have exponentially large ra- 
dial distances Ro and R from the origin. This explains the 
picture we obtain in the numerical solution and agrees with 
earlier results obtained by Eichler (1993) and by Tomimatsu 
(1994). 

In the logarithmic scale of Figs. (Pja) the impression one 
may get is that all poloidal fieldlines are focused towards 
the symmetry axis. Such a figure was produced in Saku- 
rai (1985). However, this way of plotting the shape of the 
poloidal streamlines is a deceiving about the asymptotical 
shape of the magnetosphere which is different, as seen in the 
linear scale of Fig. (^b) . In this plotting it is clear that the 
geometry of the poloidal field lines is such that a cylindri- 
cal core (the jet) is formed around the symmetry axis with 
a width of the order of the base radius, but all other field 
lines go to straight asymptotes filling all space. 



8 ACCRETION DISK-LIKE ROTATION, 

An analysis of the asymptotical behaviour of MHD outflows 
in Bogovalov (1995) shows that a discontinuity in the total 
electric current is formed on the equator for outflows from 
astrophysical objects having a magnetic field directed in op- 
posite directions in the upper and lower hemispheres. In a 
pure monopole-like magnetic field the electric current leaves 
the star in the upper hemisphere and returns back in the 
lower hemisphere. In such a magnetosphere the total elec- 
tric current is continuous on the equator. But already in the 
split monopole model where the magnetic field has opposite 
directions in the upper and lower hemispheres, the closure 
of the electric circuit occurs along the equator so that on 
the equator XB^, — 0. This means that the total electric 
current flowing in the upper or lower hemispheres is equal 
to zero. On the other hand, according to Bogovalov (1995) 
the magnitude of XB^ is not equal to zero on all field lines 
of the supersonic flow where Q ^ 0. This means that the 
function XB$(il)) shown in Fig. (|§]a) has actually a discon- 
tinuity at the point ip — 1 (on the equator) for the split 
monopole solution. This is the usual MHD contact disconti- 
nuity between the flows in the upper and lower hemispheres 
with the magnetic field equal to each other in magnitude but 
with different directions. It is important to check the validity 
of this conclusion for a pure monopole-like magnetic field. 
In the following, we shall investigate the process of the 




Figure 10. Shape of the poloidal magnetic field lines in the near 
zone of a differentially rotating magnetic rotator with a = 1: the 
field lines focus towards the pole and the eguator. 




Figure 11. Shape of the poloidal magnetic field lines for a dif- 
ferentially rotating magnetic rotator with a = 1, plotted in a 
logarithmic scale up to a radius 1000 times the initial base ra- 
dius. 



formation on the equator of such a discontinuity in the total 
electric current, in the ideal MHD approximation. For this 
purpose, we take the rotation of the star such that the total 
electric current which flows in each of the upper or lower 
hemispheres to be independently equal to zero. For exam- 
ple, this happens with a differential law of rotation wherein 
f2(V0 = on the equator and the fieldline on the equator 
does not rotate. This guarantees that XB$ — 0, everywhere 
on the equator. 
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In this work we took the simplest law of rotation 
0(V>) = a(l - ip) . 



(41) 



It is worth to pay attention to the fact that that this law 
of rotation describes qualitatively the rotation of accretion 
disks where the largest angular velocity is in the inner edge 
of the disk and the lower one at the outer edge of the disk. 
The solution of the problem for this differential law of rota- 
tion in the nearest zone is shown in Fig. (|Io|). 

The most important difference with isorotation is the 
coincidence of the fast mode and Alfvenic surfaces on the 
equator. This is simply the consequence of the assumed zero 
toroidal magnetic field on the equator. 

The poloidal magnetic field in the far zone is shown in 
a logarithmic scale in the next Fig. [ll| It is clearly seen 
in this figure that along with the expected collimation of 
the flow towards the axis of rotation, there is also a focus- 
ing of the flow towards the equator. The explanation of this 
result is simple. With the assumed differential law of rota- 
tion, the largest toroidal magnetic field in the nearest zone 
is found somewhere at the middle latitudes. The pressure of 
this toroidal magnetic field pushes the plasma towards the 
equator and the pole from those mid latitudes. 

Fig. ([jj]) shows the formation of a discontinuity in the 
total electric current near the equator (ip — 1) . On the other 
hand, iso-current contours in Fig. (fl2|b) show that at small 
r) the dependence of XB^, on ip is that of a smooth function 
with a maximum placed at ip — 0.5. This maximum moves 
to the equator as r/ increases and finally the distribution of 
XB$ is similar to the one for an isorotation as that shown 
in Fig. ^, with the exception of a discontinuity formed near 
the equator. 



9 RELATIVISTIC OUTFLOWS 

In this section we discuss briefly the problem of collimation 
of a relativistic plasma. This relativistic version of the prob- 
lem under consideration differs from the nonrelativistic one 
only in the modification of the function G(rj, ip). In the rela- 
tivistic case, the terms with the electric field play an impor- 
tant role. To specify the boundary conditions in the super 
fast magnetosonic region we will use the approximate solu- 
tion obtained in Bogovalov (1997) for a relativistic plasma 
outflow in the nearest zone. In this solution the magnetic 
flux function is 



ip = 1 — cask 



(42) 



corresponding to the poloidal magnetic field of a magnetic 
monopole and 



-xB. 



7o 
P U 



(43) 



where 70 is the initial Lorentz factor of the plasma which is 
ejected from the surface of the star, Uo is the initial four- 
velocity of the outflow, x is the cylindrical distance in units 
of the light cylinder radius. We assume a uniform rotation 
of the star, 



U P = U , 



u v =o, 



1 \ n o 



(44) 



where no is the density of the plasma on the light cylin- 
der. The solution written above is an approximate one, with 



corrections in the subfast magnetosonic region of the order 
cr/70, where a = (B 2 /4Tvmc 2 n)(R s t ar Q/c) 2 is the Poynting 
flux per particle at the equator. The fast mode surface in 
these variables is r/ = y/cr/jo- This solution is valid under 
the condition 0-/7 3 C 1. To get an idea of the physical pa- 
rameters it may be useful to recall that for the Crab pulsar 
we have a ~ 10 6 , 70 ~ 10 3 , and a/-)l ~ 10~ 3 . 

The boundary conditions in the form of the integrals 
W(ip), L(ip), Cl(ip), F(ip) were specified directly after the 
fast mode surface and the solution obtained is shown in Fig. 
(13) for sigma a = 300 and 70 = 30. The collimation of the 
flow in this relativistic case is very weak. To explain this let 
us estimate the dependance of the radius of curvature of the 
star on the distance to the star assuming that the magnetic 
field to a first approximation is the field of the magnetic 
monopole. At large distances, we can neglect the poloidal 
magnetic field and the azimuthal rotation of the plasma. In 
this case the transfield equation becomes, 



8 



8nx 2 dtp 



(x 2 (H 2 - E 2 )) 



(U p + F(j:)x 2 B p ) 
4rrxR c F(ip) 



0. 



(45) 



It is easy to obtain similarly to the nonrealtivistic case the 
approximate equation for a field line, 



dO _ sin 26 

dr ~ r[/ 2 (^ +sin 2 <9)' 

An integration of this equation gives approximately, 
~A0[/ o 2 (^ + sin 2 ( 



r f exp 



sin 26» 



(46) 



(47) 



where 77 is a lower limit of the integration. It may be seen 
that a large factor Uq/u 2> 1 is present in the expression 
above for small angles 8 < y/Uo/cr and an even larger mul- 
tiplier Uq at larger angles. Collimation is firstly expected at 
small angels A6 ~ 9. Therefore, the distance at which the 
jet is formed is 



Tcoii = Tf exp 



u± 

2a' 



(48) 



For the parameters used in the calculations r co u — 3.5 x 
10 19 r/. It is interesting to estimate the distances at which 
we should expect collimation of the wind from the Crab 
pulsar. For this pulsar a = 10 6 70 = 10 3 and thus we get 
T'coii ~ 10 226 cm, a value much larger than the size of the 
cavity (0.1 pc) formed by the wind before it terminates in 
the interstellar medium. 



This result shows that for parameters typical for radio 
pulsars there is no collimation of the supersonic plasma at 
resonable distance. Therefore, the plasma is not accelerated 
in this region as it was also concluded by Begelman & Li 
(1994). 



10 SUMMARY 

The axisymmetric 3-D MHD outflow of plasma from a mag- 
netized and rotating central object is numerically simulated 
for a wide range of angular velocities of the central star. The 
simulation of the time-dependent evolution of the flow from 
some initial state was used to obtain finally a stationary 
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Figure 12. For a differentially rotating magnetic rotator with a = 1, in (a) it is plotted the enclosed electric current XB^ by each 
poloidal field line ij) and in (b) iso-current contours in the space of (if), if). 
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Figure 13. Shape of the poloidal magnetic field lines (solid) in 
the far zone of a rotating magnetic rotator ejecting relativistic 
plasma. Dashed lines show pure radial outflow. 

solution. The model of a non rotating star with a monopole- 
like magnetic field was taken as the initial state of the mag- 
netosphere. It was found that the obtained stationary final 
state depends critically on a single parameter only. This pa- 
rameter a expresses the ratio of the corotating speed at the 
Alfven distance to the initial flow speed along the magnetic 



monopole-like fieldlines. The acceleration of the flow was 
most effective at the equatorial plane and the terminal flow 
speed depended linearly on a. Significant flow collimation 
was found in fast magnetic rotators corresponding to large 
values of a > 1, while very weak collimation occurs in slow 
magnetic rotators with small values of a < 1. Part of the 
flow around the rotation and magnetic axis is cylindrically 
collimated while the remaining equatorial part obtains ra- 
dial asymptotics. The transverse radius of the jet is found 
to be inversely proportional to a while its density grows lin- 
early with a. For a > 5 the magnitude of the speed of the 
flow in the jet remained below the fast MHD wave speed 
everywhere. We predict that a regime of nonstationary jet 
ejections may be possible at such high values of a. 

The above results have been obtained under several sim- 
plifying assumptions, such as, the neglect of gravity and 
thermal pressure as well as by taking for the initial magne- 
tosphere a split-monopole configuration. Despite those as- 
sumptions however, we recover the main results of rather 
general theoretical studies on the formation of a collimated 
outflow from a magnetized and rotating astrophysical ob- 
ject. Nevertheless, for a meaningful comparison with the 
observations one needs to relax these assumptions, a task 
taken up in the next paper. 
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